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Abstract.  Approximate  method  for  Monte-Carlo  modeling  of  interaction  of  high-altitude  engine  exhaust  plume  with 
atmosphere  has  been  developed.  The  method  was  applied  for  analysis  of  recent  experimental  data  on  Russian  spacecraft 
Soyuz  plume  glow.  In  particular,  sensitivity  of  computed  glow  radiation  to  the  choice  of  elastic  collision  cross-section  of 
atmospheric  oxygen  with  plume  molecules  and  its  reaction  cross-section  with  water  molecules  in  the  plume  is  analyzed. 
Conclusions  are  made  about  these  cross-sections  values. 


INTRODUCTION 

A  set  of  experiments  dedicated  to  measurement  of  Russian  spacecraft  (SC)  Soyuz  and  Progress  plume  glow 
in  the  region  270-340  nm  at  altitudes  of  360-400  km  was  accomplished  recently  on  board  of  Russian  Space  Station 
Mir  [1-4].  In  this  experiments  radiometric  data  on  the  intensity  and  spatial  distribution  of  the  radiation  in  the  SC 
motor  plumes  were  acquired.  A  hypothesis  was  put  forward  that  radiation  is  a  product  of  the  hypervelocity  reaction 
between  water  molecules  available  in  plume  with  atmospheric  oxygen  atoms 

0  +  H2O^OH(A)  +  OH(X)  (1) 

The  UV  glow  was  observed  in  the  region  extended  over  some  kilometers  from  the  SC  motor.  Several  papers  have 
been  published  by  now  where  this  glow  was  modeled  by  DSMC  method  [5-7].  Results  of  the  modeling  overpredict 
the  experimental  data  on  one  or  two  orders  of  magnitude,  depending  on  the  model  parameters.  Conclusion  was  made 
that  main  source  of  this  discrepancy  might  be  not  appropriately  chosen  model  parameters,  such  as  elastic  and 
reaction  (1)  cross-sections.  No  direct  measurements  or  accurate  computation  of  the  cross-section  values  in  the  range 
of  collision  energies  relevant  to  the  experiment  are  known  so  far.  So,  for  modeling  of  these  cross-sections  authors 
used  VHS  and  TCE  models  correspondingly  and  the  model  parameters  were  derived  by  extrapolating  of  the 
parameters  defined  at  low  energetic  collision  to  the  high  energetic  collision  case.  VHS  model  parameters  were  used 
as  those  obtained  from  the  viscosity  measurements.  Whereas,  parameters  of  TCE  model  were  derived  from 
measured  rate  constant  of  the  reaction  (1),  but  resulting  in  not  excited  products 

0  +  H20^20H(X) 

Such  approach  was  not  completely  satisfying  from  the  standpoint  of  the  comparison  with  experimental  data. 
Therefore,  necessity  arises  to  analyze  the  glow  phenomena  in  more  detail  and  revise  more  accurate  collision  cross- 
section  values  basing  on  the  newly  obtained  experimental  data.  This  analysis  requires  a  vast  amount  of  parametric 
computation  and  can  not  be  effectively  accomplished  basing  on  the  time  consuming  full  scale  DSMC.  As  an 
alternative,  a  simplified  physically  transparent  model  of  the  glow  was  developed  for  the  data  analysis.  In  this  model 
averaged  elastic  cross-section  values  are  used  to  compute  the  flow-field  and  overlay  method  is  employed  for 
radiation  computation.  In  comparison  with  full  DSMC  the  simplified  computation  requires  much  less  time  resulting 
in  the  same  outputs.  First  attempts  to  apply  this  model  to  analysis  of  the  experimental  data  fully  justified  the 
approach. 
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2.  PHYSICAL  AND  COMPUTATIONAL  MODEL 


2.1.  Input  data  for  model  statement. 

The  mathematical  model  is  intended  for  the  analysis  of  concrete  experimental  data.  Therefore  standard 
conditions  of  realization  of  experiment  are  presented  below.  The  common  scheme  of  observations  is  shown  on  Fig. 
1 .  Two  schemes  of  observation  were  applied:  a)  observation  «  in  a  tail  »,  at  a  small  viewing  angle  /J ,  and 
observation  from  the  side  at  tangent  angle  /3  =  90°.  In  both  schemes  the  pitch  angle  (X  did  not  exceed  1 7°. 


FIGURE  1.  Scheme  of  space  experiment. 

The  MIR-SOYUZ  range  varies  from  1 1  km  up  to  25  km.  For  numerical  simulation  the  standard  values  of  parameters 
of  a  trajectory  and  atmosphere  were  selected: 

77=360-400 km,  F=7.35km/s,  na  =  1.3410l4m~3  ,  7=886 K, 

X(O)=0. 909,  X(He)=0.05n,  X(N2)  =  0.0393  . 

Here  X(M)  is  a  mole  fraction  of  M  gas  component. 

As  the  initial  data  for  calculation  of  rarefied  flow  of  gases  the  data  obtained  at  solution  of  continuum  equations 
(up  to  distance  about  20  meters  from  the  nozzle  exit  section)  were  used.  The  following  values  of  parameters  were 
adopted  on  plume  axis  at  x  =  20.18 m:  density  n  =  0.331  10  m  ,  gas  velocity  (in  a  nozzle  coordinate  system) 
u  =  3335 mis  and  temperature  T  =  37.52/1 . 

It  was  supposed,  that  the  signal  of  UF-imager  mainly  was  defined  by  the  contribution  of  an  emission 
OH(A  —>  X) .  The  characteristics  of  reaction  (1)  are  those:  threshold  energy/^  =  4.S3eV,  lifetime  of  excited 

molecules  T  =  1 0  9  s  and  average  energy  of  quantum  hv  -  6.43  10  19  J. 

2.2.  Mathematical  model 


2.2.1.  Simple  physical  model. 

The  input  data  and  notion  about  emission  source  allow  to  present  the  simplified  physical  model.  Let's  consider 
plume  in  a  spacecraft  coordinate  system,  with  the  nozzle  as  an  origin.  It  is  well  known,  that  density  in  plume  is 

proportional  to  (x/xs)~2 ,  where  x  -  distance  up  to  a  nozzle.  This  dependence  enables  to  estimate  the  distance,  on 
which  numerical  density  in  plume  is  compared  to  numerical  density  of  gas  of  a  nonperturbed  atmosphere: 
Lx  =10  km  .  Such  estimate  displays  the  weak  influence  of  an  atmosphere  on  a  plume  field  at  a  distance  of  about  1 
km.  On  distances  comparable  with  Lx ,  it  is  necessary  to  take  into  account  the  interaction  of  a  plume  and  atmosphere. 

Temperature  of  a  plume  is  so  low,  that  it  is  possible  to  neglect  by  thermal  velocity  of  plume  molecules.  The 
threshold  of  process  (1)  is  high.  Therefore  it  happens,  in  main,  owing  to  existence  of  oxygen  atoms  thermal  velocity. 
If  the  oxygen  atom  has  collision  with  a  plume  molecule,  then,  first,  it  removes  from  area,  close  to  plume  axis,  and, 
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secondly,  its  velocity  decreases.  Therefore  the  energy  of  relative  motion  of  oxygen  atoms  and  molecules  of  water 
becomes  lower  energy  of  a  threshold.  Therefore  at  x  <  Lx  it  is  admitted  to  consider  following  model  of  gas  flow. 
The  creation  of  a  plume  at  x  «  L{  happens  irrespectively  to  interaction  with  an  atmosphere.  The  atoms  of  oxygen 
scattered  on  plume  molecules.  As  the  thermal  velocity  of  oxygen  atoms  is  much  less  than  their  average  velocity,  it  is 
possible  to  accept  the  collision  cross  section  Oe  to  be  constant  in  the  flow.  The  generation  of  OH  {A)  results  from 
O  +  H20  collision  according  to  probability  of  reaction  (1). 

2.2.2.  Mathematical  model. 

The  inelastic  collisions  frequency  of  water  molecules  with  oxygen  atoms  which  leads  to  emission  is  equal 

W(x)  =  JJ  f„  (f „)/,(! x)gor  {E)dladl „  E>Eth,  (2) 

Here  f„,f\  are  the  velocity  distribution  functions  of  oxygen  atoms  and  water  molecules,  or(E)  -  inelastic  collision 
cross  section.  In  accordance  to  2.2.1.  these  functions  are  the  solutions  of  simplified  Boltzmann  equations: 

z*P-=yL-j(fi’f]'>  (3) 

dxk  '-‘J 

$k^  =  -foo^figai(.L4i)4i~-f(4  o-h'V  (4) 

Here  /00  is  the  distribution  function  of  oxygen  atoms  which  move  without  collision  with  plume  molecules,  (Je  - 
effective  elastic  cross-section,  and  it  was  taken  into  account,  that  the  change  of  relative  concentrations  of 
components  in  plume  may  be  neglected.  The  equation  (3)  was  solved  by  the  DSMC  method  and  linear  equation  (4) 
by  the  usual  Monte-Carlo  method.  The  solution  of  Eq.(3)  have  shown  that  the  density  of  plume  has  the  similarity. 

2.2.3.  Collision  cross-section  adopted  values. 

Elastic  cross-section.  The  data  on  intermolecular  potentials  of  oxygen  atoms  and  plume  molecules  [8-10]  are 
sparse.  The  most  of  data  were  obtained  in  experiments  with  scattering  beam  [8]  for  all  components  of  a  jet  excepting 
H20 .  Let  us  compare  the  values  of  transport  collision  cross-section  Q  for  O-  N2  collision,  calculated  with  the 

data  of  [8,10],  At  average  velocity  of  oxygen  atoms  these  values  are  7.42  10  20  m 2  and  7.96-10  20  m2  accordingly. 
The  estimate  uses  the  Bom-Mayer  exponential  repulsive  potential.  The  cross-sections  calculated  with  the  help  of 
Lennard-Jones  potential  are  approximately  twice  more.  The  similar  results  may  be  obtained  for  the  O  -  CO  and 
0-C02  collision  cross-section.  These  estimates  lead  to  parametric  investigation  at  cross-section  value  range 

7.5-10  20 m2  <  ae  <  20-10  20 m 2 .  Nevertheless  authors  suppose  by  more  natural  in  the  considered  problem  to  use 
the  cross-section  values  obtained  in  beam  experiments. 

Inelastic  cross-section.  The  data  on  the  process,  interesting  for  us,  practically  miss.  The  collision  model  was  used 
conventional  for  simulation  of  inelastic  processes  [11],  though  this  form  may  be  inadequate  to  the  true  process: 

°  r=°  0i(E~ EthJ)n‘  E~\e>  EthJ  (5) 

_ 1/9 

Thus  two  sets  of  parameters  values  were  used:  nx  =0.5,  E^j  =4. 83eK,  <701  =1.15-10  m  (eV)  ; 
«i  =1.8,  E  th  |  =  4.79eK,  cr0i  =3.1-10 _20  m 2  (e  V) -0'8 .  The  relation  (5)  closes  a  statement  of  mathematical  model  of 

processes  (2),  (3),  (4),  which  leads  to  emission.  This  model  enters  splitting  processes  and  allows  to  calculate 
scattering  and  inelastic  collisions  of  oxygen  atoms  after  the  parameters  of  a  plume  are  defined. 

2.2.4.  Analytical  model. 

Interaction  of  oxygen  atoms  with  plum  molecules  (see  Eq.  (4))  in  case  of  a  cylindrical  symmetry  of  may  be 
described  especially  simply  on  axis.  In  this  case  it  is  possible  approximately  to  split  the  description  of  processes 
leading  to  emission.  Namely,  to  divide  the  process  of  elastic  scattering  of  oxygen  atoms  and  their  inelastic 
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interaction  with  molecules  of  water.  Let's  assume,  that / 00(x,<f;)  =  «(x)/  oo,«(£)  ■  Integrating  Eq.(4),  substituting 
outcome  in  Eq.(2)  and  using  a  smallness  of  SZl ,  we  can  receive: 

W(x, r  =  0)  =  /ae )X(H20)Xm  ( 0)naV  •  x;lKnplfi  ■  y2  exp(-y)  (6) 

y  =  Kn'pifi ix*  /*)  >  Kn'l 0  =(V  +  u)  V~1aeA-  nplfi  ■  x* , 

°r,eff  =  (*-  (O)  ■  na  r1  (V  +  «)_1Jar  (£)g  •  /eo  (oo)rff  (7) 

So  simple  computational  model  allows  to  make  the  following  conclusions.  The  maximum  value  of  inelastic 
collision  frequency  is  reached  at  y  =  2,  i.e.  at  point  xmax ,  distance  from  which  to  nozzle  is  proportional  ae . 
Maximum  value  is  in  inverse  proportion  to  oe .  If  we  observe  the  emission  from  infinitely  large  distance  along 
a  centerline  the  radiance  will  be  peer 

Wx  =  ^  f  W(x^dx  =  T~  X pi  (^20)  •  iP)  n  V ^  (8) 

4  k  Jo  4  n  ae 

It  is  seen  that  the  radiance  is  in  inverse  proportion  to  oe . 

2.3.  The  model  defects  assessment. 


2.3.1.  Error  from  the  limitation  of  a  computational  area. 

For  the  numerical  modeling  of  the  process  (4)  the  conditions  of  a  undisturbed  flow  are  put  on  final  distance  from 
a  nozzle  Lmm  .  In  this  case  in  expression  (8)  there  is  a  factor  krja :  =  exp(k'n“/10(xj,/Z,max)).  The  value  of  this 

coefficient  at  x*  =19.18  m,  Lmax  =15  km,  ae  =15TO“20m2  is  equal  to  1.183.  Thus,  the  limitation  of  an  axial  size 
of  computational  area  by  15  km  reduces  the  outcome  near  to  an  axis  up  to  18  %. 

2.3.2.  Error  caused  by  neglecting  of  plume  molecule  scattering. 


Consider  the  process  of  plume  molecule  scatter  at  large  distance  from  a  nozzle  in  the  same  way,  as  the  oxygen 
atoms  scattering.  It  is  easy  to  receive  a  following  set  of  equations  for  the  plume  numerical  density  on  an  axis  of  a 
plume: 


ix  npd  =  -n pin a{V  +  u)oe,  V  —  na  =nplna{V  +  u)oe 
x  ox  d  x 


(9) 


The  solution  of  equations  (9)  by  the  method  of  matched  asymptotic  expansions  on  the  small  parameter, 
(eKnpl  0) , £  =  (V  +  u)u^[ x.ro e  ana  ,  we  shall  receive  in  main  approximation:  Wl/W2  =  exp(-ex) .  Here  W1  -  the 


value  of  radiance  on  axes  appropriate  to  model,  in  which  it  is  taken  into  account  the  plume  scatter,  W2  -  the  value 
appropriate  to  model  without  the  scatter  account.  The  ratio  Wxi  /  Wx2  under  the  observation  from  distance  of  15  km  is 
equal  0.877  and  0.785  for  oe  a  =  10  10  20  m1  and  1510_20w2  accordingly.  Thus,  the  neglecting  by  the  plume 

scattering  leads  to  an  error  about  20  %.  Let's  remark  also,  that  the  errors  of  models  caused  by  limitation  of 
computational  area  and  neglecting  by  the  plume  scattering  have  different  signs  and  partially  cancel  each  other. 


SIMULATION  RESULTS. 


3.1.  Models  comparison. 

The  equations  (2)-(5)  of  mathematical  model,  were  solved  by  a  Monte-Carlo  method.  The  results  obtained  with 
the  help  of  this  model  (curve  1)  are  compared  in  a  Fig.  2a  to  the  analytical  profile  JV(x,0)  (see  eq.  (6)  -  curve  2)  and 
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FIGURE  2.  Inelastic  collision  frequency  profiles:  a)  comparison  of  computational  models,  b)  dependence  on  elastic  collision 
cross-section. 

with  complete  DSMC  result  [6]  (points  3).  The  normalized  form  Wnorm  =  W /  Wmax  and  xnonn  =  x/ xmax  is  chosen 
for  comparison  in  accordance  with  Eq.(6)  since  it  should  not  depend  on  uncertainties  in  choice  of  values  ae,ar  eff  . 
The  good  agreement  of  the  profiles  testifies  that  the  adopted  simplifications  of  model  are  justified. 

3.2.  Analysis. 


3.2.2.  Dependence  on  cross-section. 

The  analytical  model  not  only  correctly  determines  the  form  of  the  profile,  but  also  predicts  dependencies  of  the 
values  and  xmax  from  the  value  of  elastic  collision  cross-section.  In  a  Fig.  2b  the  profiles  of  W  on  an  axis  of 

plume  presented  for  axially  symmetric  flow  for  three  values  of  collision  cross-section  ae  =ae0  \  0  20  m2 , 
l-cre0  —10,  2  —  oe0  =15,  3-CgQ  =  20 ,  and  inelastic  cross-section  (7H  (see  2.2.3.).  It  is  easy  to  test,  that  the  detected 
earlier  approximate  dependencies  for  a  position  of  maximum  xmax  and  its  value  Wmdx  from  value  of  ae  are  fulfilled. 

In  Fig.  3a  the  profiles  of  radiance  on  an  axis  of  a  plume  are  shown  at  an  observation  from  the  side 
perpendicularly  axes  of  plume,  W,  (x,y  =  0)  =  (hv/ 4/r)J  Wdz ,  for  the  same  three  values  of  elastic  collision  cross- 

section.  From  this  data  it  follows,  that  the  effective  value  of  elastic  collision  cross-section  between  atmospheric 
oxygen  atoms  and  plume  molecules  may  be  detected  by  results  of  glow  observation  from  the  side  provided  that  the 
inelastic  processes  do  not  influence  the  form  of  measured  profiles.  The  absence  of  such  influence  is  illustrated  in 


(a) 


FIGURE  3.  Radiance  profiles:  a)  side  view  Wz(x,y=0)  -  dependence  on  elastic  collision  cross-section,  b)  tail  view  WS(R)  - 
dependence  on  inelastic  collision  cross-section. 
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Fig.  3b:  the  relative  profiles  of  the  value  WX(R)  (R-distance  from  axis)  for  two  cross-section  of  the  inelastic 

process,  see  Eq.  (5),  are  shown  at  ae  =15- 10  20  m2  .  Difference  in  the  profiles  is  insignificant.  These  data  justify 
the  splitting  processes  on  two,  carried  out  by  analytical  model  (scattering  of  fast  atmospheric  oxygen  atoms  and 
inelastic  collisions  of  atoms  with  molecules  of  water)  not  only  on  an  axis  of  plume. 


FIGURE  4.  Dependence  of  Wx  on  elastic  collision  cross-section. 

In  Fig.  4  the  profiles  of  the  value  Wx ,  i.e.  dependence  of  this  value  on  distance  from  an  axis,  are  shown  at  the 
same  three  values  of  elastic  cross-section.  The  values  of  collision  cross-section  influence  the  profile  of  this  value, 
but  only  at  small  distance  from  an  axis.  At  increase  of  distance  this  influence  disappears.  Such  dependence  is 
explained  by  the  fact,  that  outside  of  plume  core  the  flux  of  oxygen  atoms  does  not  damp. 

3.2.3.  Comparison  with  experiment. 

These  results  give  the  possibility  to  define  the  value  of  effective  collision  cross-section  cre  not  only  by  the 
measurements  of  Wz  but  also  by  the  measurements  transverse  profiles,  the  value  Wx .  The  result  of  such 
determination  of  the  ae  value  is  shown  in  Fig.  5.  The  normalized  profiles  WLnorm  observed  from  the  finite  distance 
(see  Fig.  1)  were  used  for  this  purpose,  because  it  does  not  depend  on  or .  The  angle  x  is  the  angle  between 
observation  line  and  MIR-SOYUZ  line.  Curves  1  and  2  correspond  to  the  computation  with  rre=7.5-10'20  m2  and 
ae=  15-1 0  20  m2.  Profile  derived  from  the  experiment  is  also  shown  in  the  figure.  The  experimental  data  are  evidently 
closer  to  curve  1,  that  confirms  oe  estimate  obtained  from  beam  scattering  measurement  data  (see  2.2.3.)  At  this 
computation  observation  geometry  of  real  experiment  performed  on  28  February  1999  was  used  and  non-symmetry 
of  the  flow  was  taken  into  account  as  well. 


FIGURE  5.  Comparison  of  experimental  and  numerical  normalized  radiance  profiles. 
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Figure  6  gives  an  example  of  numerical  modeling  of  radiation  map  observed  in  the  same  experiment.  This 
modeling  revealed  that  maximal  brightness  in  the  map  changed  not  very  much  while  the  observation  condition  and 
the  flowfield  parameters  were  varied.  Assuming  the  atmosphere  parameters  and  the  measured  radiation  spectral 
content  are  known,  one  could  derive  the  oreff  value  from  these  maps.  The  first  result  of  so  comparing  of 

experimental  and  numerical  data  leads  to  the  upper  limit  of  effective  inelastic  collision  cross-section 

—22  2 

o  reff  =0.9  10  m  .  This  value  close  to  the  water  photo-dissociation  cross-section  with  OH  {A)  coming  out  as 
product.  The  received  value  <7 r  eg  is  lower  than  the  values,  computed  by  Eq.(5),  by  factor  of  10-25. 


These  examples  shows  that  suggested  method  can  be  effectively  used  for  analysis  of  the  flowfield  and  derivation 
of  the  cross-section  values  from  the  experimental  data.  It  should  be  noted  that  accuracy  of  these  derivation  may  be 
affected  by  not  only  the  experimental  errors,  but  by  adequacy  of  the  model  approximation  as  well.  For  instance, 
or  eg-  derived  from  the  experimental  data  evidently  depends  on  the  degree  of  water  condensation  in  the  rarefied 
plume.  This  fact  was  ignored  in  the  frame  of  the  model.  In  future  we  hope  to  improve  physical  model  of  the 
phenomena  and  obtain  more  accurate  estimate  of  the  <3r,eff  value.  And  finally  -  the  method  suggested  here  may  be 

used  for  derivation  of  the  energy  dependence  of  ar  (E)  by  solution  of  inverse  problem  (see  2.2.4.).  This  solution  will 
require  incorporation  of  various  ar  ^  measured  at  different  angles  between  the  plume  axis  and  atomic  oxygen  flux 
velocity. 
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